Direct Estimate of the Static Length-Scale Accompanying the Glass Transition 
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Characterizing the glass state remains elusive since its distinction from a liquid state is not obvious. 
Glasses are liquids whose viscosity has increased so much that they cannot flow. Accordingly there 
have been many attempts to define a static length-scale associated with the dramatic slowing down 
of supercooled liquid with decreasing temperature. Here we present a simple method to extract the 
desired length-scale which is highly accessible both for experiments and for numerical simulations. 
The fundamental idea is that low lying vibrational frequencies come in two types, those related to 
elastic response and those determined by plastic instabilities. The minimal observed frequency is 
determined by one or the other, crossing at a typical length-scale which is growing with the approach 
of the glass transition. This length-scale characterizes the correlated disorder in the system, where 
on longer length-scales the details of the disorder become irrelevant, dominated by the Debye model 
of elastic modes. 



The phenomenon of the glass transition in which super- 
cooled liquids exhibit a dramatic slowdown in their dy- 
namics upon cooling is attracting a tremendous amount 
of effort. One of the thorny issues has to do with at- 
tempts to define and measure a static length-scale that 
will grow upon the approach to the glass transition. Sev- 
eral such attempts were published in the recent litera- 
ture, using delicate measures of higher-order correlation 
functions T, '5] , of the effects of boundary conditions, 
point to set correlations Q, the scaling of the non afhne 
displacement field Q and patch correlation scale @. In 
this Letter we propose a different approach which yields 
a very natural static length that appears to fit the bill. 
An advantage of this approach is that it can be imme- 
diately and easily applied to both experimental fio'l and 
numerical data (see below). 

Our starting point is the fact that at lower frequency 
the density of state (DOS) for a disorder system refiects 
the access of plastic modes which are not exhibited by 
the de nsity of states of the corresponding crystalline ma- 
terials This excess of modes is sometime referred 
to as the 'Boson Peak' It is most natural to discuss 
the 'modes' in terms of the eigenfunctions of the Hessian 
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where U{ri,r2, ■ ■ ■ rjv) is the total energy of the system 
consisting of TV particles, and {rk}^^i are their coordi- 
nates. Being real and symmetric, the Hessian is diagonal- 
izable, and in all discussions below we refer to this matrix 
with the Goldstone modes (zero modes due to symme- 
tries) being pruned out. Recently we had progress in 
understanding the analytic form of the density of states 
of generic glassy systems IJj. The eigenvalues appear 
in two distinct families, one corresponding to the Debye 
model of an elastic body, while the density of access plas- 
tic modes could be approximated as 

/pl(A) = B(T)A^for A^O, (2) 

where the pre- factor B{T) being strongly dependent on 



temperature and the exponent 9 being weakly dependent 
on temperature. This dependence is a partial measure 
of the degree of disorder which grows with temperature, 
or for quenched systems, is partially controlled by the 
quench protocol from the melt to the glass. Together 
with the standard Debye contribution one can approxi- 
mate the low-frequency tail of the density of states with 
the following approximation 
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The first term , which is the Debye approximation of 
the DOS for a continuum elastic medium, contains Xd — 
^p2/d-i ^ which is the Debye frequency and /i is the shear 
modulus. For the sake of generality we do not assume 
here any form for the distribution /pi, but we will see 
bellow that the form ^ will reappear naturally in the 
analysis below. 

The main idea which will fiash out the static typical 
scale is that the minimal eigenvalue Amin observed in a 
system of TV particles will be determined by the first or 
the second term in Eq. [3l There will be a value of N 
where a cross-over will occur, such that for larger system 
the Debye form determines Amin whereas for small N the 
plastic density of states determines Amin- This cross-over 
will be interpreted below in terms of a typical length-scale 
separating correlated disorder from asymptotic elasticity. 

To see the cross over integrate Eq. ([3]) from zero to 
(Amin)- Since every quenched system is random, every 
realization will have a somewhat different value of Amin, 
and we therefore consider an ensemble average over many 
realization of system of the same number of particles N 
quenched with the same protocol. 
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FIG. 1: Left Panel : Minimal eigenvalue calculated for a model glass former interacting via a pure repulsive pair-wise potential 
plotted as a function of system size. Right Panel: Minimal eigenvalue rescaled by the high frequency shear modulus ( which is 
proportional to the square of the Debye frequency) plotted as a function of system size. Notice that at larger system size the 
rescaled eigenvalue becomes independent of temperature indicating the crossover to the continuum elasticity limit. 



Introducing Eq. ([3]) into the integral we find 
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we end up with 
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Changing slightly this equation, writing N — pV, and 
renormalizing the coefficient A accordingly, we write 
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Finally, since the function Q is monotonically increasing, 
) can be written as a scaling function of the form 
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where T = G'^ &nd ^'^{T) 



pB(T)- 



The typical scale 



5(r) will be calculated by demanding that all the data 
calculated for different system sizes and temperatures 
should collapse into a master curve just by appropriately 
choosing the ^ (T) . 

The analysis is then done as follows. Two typical 
glass formers were considered, both obtained from a bi- 
nary mixture, one with a purely repulsive potential and 



the other with repulsive and attractive contribution (the 
Kob- Andersen model). The details of the potentials can 
be found in 0, Hj . The systems were equilibrated at some 
temperature T > Tg and then instantly quenched by di- 
rect energy minimization to the nearest local minimum 
of the energy landscape. At this state the Hessian was 
computed and the minimal eigenvalue was obtained using 
the Lanczos algorithm [oj . For a given system size N and 
temperature T this procedure was repeated to have an 
average (Amin) until convergence was achieved. At this 
point the temperature or the system size were changed 
and the procedure was repeated, to eventually have a ta- 
ble of (Amin)(A^, T). The density was fixed for every type 
of potential, (0.85 for pure repulsion in 2D, 0.82 for pure 
repulsion in 3-D, and 1.20 for the Kob- Andersen model). 

In the left Panel of Fig[Tl we have plotted the minimal 
eigenvalue (Amin) of the two-dimensional purely repulsive 
super-cooled system as a function of system size N for dif- 
ferent temperatures. In the right Panel the same minimal 
eigenvalue was rescaled by the characteristic Debye value 
Ad and plotted as a function of system size. Note that at 
larger system size they become independent of temper- 
ature consistent with the fact that at larger system size 
one recovers trivial scaling predicted by the continuum 
elasticity theory. 

In Fig. [21 we re-plotted the same data of Fig. [T] accord- 
ing to the scaling ansatz Eq. [HI The left panel exhibits 
the results without rescaling by £,'^{T). In the right panel 
we extracted the disorder length scale ^ (T) by collapsing 
the data. The resulting data collapse into a straight line 
indicates that Eq. 1^ is obeyed to high precision. The 
collapse itself supports the scaling ansatz Eq. [SI In the 
inset at the right panel we show how the typical scale in- 
creases when the glass transition is approached. In Fig,[3| 
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FIG. 2: Scaling of the minimal eigenvalue calculated for a 
model glass former interacting via a pure repulsive pair-wise 
potential in two dimensions, here and in all the figures below ^ 
was chosen as ^ = = 1 at the highest available temperature. 



FIG. 5: The scaling of configurational entropy with the length 
scale extracted from the finite size scaling of minimal eigen- 
value for the Kob- Anderson Model. The data is taken from 
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FIG. 3: The measured minimal eigenvalue as a function of 
system size for a binary system interacting via a pure repul- 
sive potential in 3 dimension. Every data point represents an 
average over 2000 Inherent structures. 



and 13] we show the same the same analysis for the purely 
repulsive glass in 3D, a and for the Kob- Andersen model 
in 3D. Both the data collapse and the resulting increase 
in the typical scale appear very encouraging and support 
the approach proposed above. Note that the approx- 
imate linearity in the collapsed scaling function in the 
right panel of Fig. [3] again indicates the relative accu- 
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FIG. 4: The measured minimal eigenvalue as a function of 
system size for a binary system interacting via Kob- Andersen 
potential in 3 dimension. Every data point represents an av- 
erage over 2000 Inherent structures. 



racy of the functional form ([5]) for the density of plastic 
modes. We have written the value of appropriate expo- 
nent 9 in Eq. ([2]) directly into the graphics. In the right 
panel of Fig. |4]we see some curvature that indicates that 
higher order terms in the density of states already play 
some role. Nevertheless the data collapse is superb. 

Finally, it is interesting to ask how our typical scale 
helps in understanding other measures of disorder that 
were proposed in the past. As an example we con- 
sider here the configurational entropy Sc{T) of the Kob- 
Andersen model in 3D. This entropy was computed in 
Refs. 0, [3] and the reader is referred to these publica- 
tions for a full description of the method and the results. 
The configurational entropy is expected to characterize 
the degree of disorder in the supercooled liquid. Once 
Sc{T) was computed, it was discovered phenomenolog- 
ically that the data at different temperatures could be 
collapsed by rescaling the system size by a typical scale 
({T). It remains however mysterious what that length- 
scale might be, and how to measure it independently of 
the configurational entropy. At this point we can offer 
a resolution of that mystery. In Fig. [5] we show, in the 
left panel, the configurational entropy measured at differ- 
ent temperatures as a function of the system size. Since 
the number of accessible minima reduces as a function 
of temperature, (high energy 'states' are excluded when 
the temperature decreases), the configurational entropy 
goes down as seen in the left panel of Fig. \5\ In the right 
panel we exhibit Sc{T)/ Sc{N — ^ oo) as a function of the 
re-scaled system size N/£^^{T). The rescaling is achieved 
using our data for ^(T) as shown in the right panel of 
Fig. |4] without any re-fitting. The data collapse shows 
that the mysterious £{T) is precisely our typical scale as 
propose in this Letter. 

We should stress that the estimate of our typical scale 
is easily achieved also in appropriate experiments. It was 
shown in Ref. that the density of states can be calcu- 
lated in experiments on colloidal glassy systems, and in 
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FIG. 6: The logarithm of the relaxation time Ta{T) measured 
for all the three models discussed above, plotted against ^'^/T. 
The straight lines are a guide to the eye. 



particular Amin is available. We propose that such mea- 
surements should be repeated as a function of system size 
at different values of the packing fraction. Such data can 
be used as explained above to determined the dependence 
of our typical scale on the packing fraction, throwing new 
light on the interesting physics of these complex systems. 

Finally, it is impossible to end this Letter without ask- 
ing how the newly found length-scale correlates with the 
dramatic slowing down that is observed in our super- 
cooled liquids as the temperature decreases. If we take 
the point of view that the typical size of relaxation events 
also depends on ^, either proportional to ^ if they are 
stringy in nature or on if they are planar, then we 
expect the relaxation time Tq, to be of the order of 



Ta OC exp 



where C is some unknown scale and t/; is 1 or 2 depending 
on the geometry of the relaxation events. We see that Eq. 
(1101) is validated with ip = I for the repulsive model in 
both 2D and 3D, and with ip = 2 ioi the Kob- Andersen 
model. It is not known at this point why indeed "0 niay 
differ when the attractive part of the potential is added, 
and we must leave this issue for further research in the 



future, and see Ref. [16|, |17[ for some comments on this 
issue. 

In summary, we have proposed here a very simple 
method to extract the typical scale that separates a dis- 
order dominated regime from an elastic dominated be- 
havior. All that is needed is the measurement of the 
minimal eigenvalue of the Hessian matrix (or, equiva- 
lently, the minimal harmonic frequency of the system), 
for systems of different size and temperature. A simple 



re-plotting procedure of the data is then used to extract 
the typical scale. This scale appears to properly collapse 
the data of the configurational entropy, resolving a rid- 
dle that existed for some time regarding the nature of the 
length-scale that does it. Finally, and only in a manner of 
passing, we also considered the relation of the obtained 
length-scale to the relaxation time of the super-cooled 
liquids, and found a strong indication that this length 
scale also determines the observed dynamics. We trust 
that the newly proposed length-scale would be computed 
in further numerical and laboratory experiments by other 
groups to enhance the understanding of the glass tran- 
sition. In particular it would be useful to compare this 
length-scale to other length-scales that were proposed by 
other groups as mentioned in the introduction. 
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